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The 3+1 formulation of scalar-tensor theories of gravity (STT) is obtained in the physical (Jordan) 
frame departing from the 4+0 covariant field equations. Contrary to the common belief (folklore), 
the new system of ADM-like equations shows that the Cauchy problem of STT is well formulated (in 
the sense that the whole system of evolution equations is of first order in the time-derivative). This 
(3JT)' is the first step towards a full first order (in time and space) formulation from which a subsequent 

~*~t . hyperbolicity analysis (a well-posedness determination) can be performed. Several gauge (lapse and 

. shift) conditions are considered and implemented for STT. In particular, a generalization of the 

harmonic gauge for STT allows us to prove the well posedness of the STT using a second order 
analysis which is very similar to the one used in general relativity. Some spacetimes of astrophysical 
and cosmological interest are considered as specific applications. Several appendices complement 
■ the ideas of the main part of the paper. 
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In the last decade there has been an increasing interest in the so-called scalar-tensor theories of gravity (STT; see 
Ref. [lj for a review) in view of the possible deviations that Einstein's general theory of relativity (GR) could show 
£h ■ in the framework of several upcoming observations. The most prominent example of STT is perhaps the Brans-Dicke 
77' (JD) theory which intends to incorporate the Mach's principle into GR by considering a varying gravitational constant 

>: 0. 

Scalar-tensor theories are alternative theories of gravity in which a scalar-field is coupled non-minimally (NMC) to 
the curvature (in the Jordan-frame representation -see Sec. II-). Due to this fact and also due to the mere presence 
of a fundamental scalar field, these kind of theories predict several phenomena not present in GR. In particular, there 
has been a considerable interest in detecting a scalar-wave component (spin-0 waves) in addition to the ordinary 
gravitational waves (spin-2 waves) predicted by GR (see for instance Refs. yl 0,0)- Among the potential emitters of 
such scalar-waves are the compact binary systems and neutron stars. On the other hand, if no appreciable deviation 
of GR are observed, then the experiments will be equally helpful to bound the STT parameters or couplings In, B , Bl ■ 
Through the effects of spontaneous scalarization in neutron stars (a new phenomenon predicted by STT) |6L l9t. 1 1 OL 1 1 1| . 
it seems that scalar-waves have a chance to be observed in detectors like VIRGO or LIGO within a few hundreds 
of kiloparsecs [13 . Moreover, there is the hope that resonant mass detectors of spherical shape or interferometers 
like LISA can resolve the existence of a scalar-gravitational wave 0,[3- Furthermore, scalar- waves produced during 
oscillations in neutron stars could leave an imprint on the spin-2 waves spectrum and therefore spin-0 waves could be 
indirectly detected |12|, 

At the large scale (cosmol ogic al scales), STT have been proposed as models for dark energy that can replace the 
cosmological constant 0, 0, [HJ , and also as an attempt for explaining some other features of the galaxy distribution 
in our universe [l6|. 

At smaller scales (scales of the order of meters or kilometers) the spurious discovery of a fifth force ^tJj renewed 
the interest about the existence of new fundamental fields of meter-range ^3- At this regard, considerable effort has 
been put in measuring gravitationally such kind of interactions that can mimic a varying gravitational "constant" G. 
Indeed, a short range scalar field in STT would induce deviations of the 1/r 2 Newtonian law due to the appearance 
of a Yukawa type of interaction (see Appendix B). For instance, at sub-millimeter scales, experiments using torsion- 
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balance instruments bound such a Yukawa type of gravitational force |l9j. Other tests consist in the analysis of 
gravitational signals induced by variations on the mass of a lake |20j| . Such experiments probe basically fields within 
the range of meters to some kilometers. Another kind of experiments probes variations of the Newton's gravitation law 
at scales of two Earth's radii by measuring the gravitational effects on the orbit of the laser-ranged LAGEOS satellite 
[2lj| . At solar-system scales, several experiments bound the deviations of STT from GR through the post- Newtonian 
parameters (see Ref. [2^). In particular, the Cassini probe restricts the effective BD parameter wbd > 4 x 10 4 pjjj . 

At scales ~ 10 Mpc, large-scale structure surveys constrain also the strength of a Yukawa type of gravitational 
interaction |24j . 

The WMAP data [2j|, can in fact be used to bound the variation of the gravitational constant in some class of 
STT 0. 

In view of the above research, STT are very well motivated theories from the phcnomenological point of view. On 
the other hand, it is well known that more fundamental theories like strings, M-theory, etc., predict in the low energy 
limit the existence of scalar-fields that can be coupled non- minimally to the curvature (see Ref. [2(| for a review) . 

From the mathematical point of view it is important to provide a solid basis for STT. For instance, a meaningful 
physical theory is expected to have a well formulated and well posed Cauchy problem (initial value problem) . A 
formidable attention to this issue has been given to pure GR, from the pioneer investigations of several mathematicians 
(e.g. see Refs. [2^) until the present where the task of analyzing hyperbolic formulations of GR in several gauges 
have been renewed (see Ref. 28] for a review) . Hyperbolic systems are important for numerical stability. This in turn 
is relevant for the numerical analysis of astrophysical systems (e.g., gravitational core collapse, black hole collisions, 
inspirilling binaries) whose effects are to be confronted with the observations. 

In the case of STT, the Cauchy problem has not been intensively analyzed like GR. The reason is that STT when 
expressed in the so-called Einstein frame (cf. Sec. II), give rise to field equations which take a very similar form as 
in GR 1], and therefore the well-posedness of the Einstein's field equations (in a given gaug e) could in principle be 
translated to STT. In the case of very particular STT like the one proposed by O'Hanlon |29f. it has been shown that 
they can be mapped to higher order theories of gravity where the Cauchy problem was studied [23] • In the Jordan 
frame, the well posedness of the Cauchy problem has been analyzed for two particular examples of STT. One is the 
BD theory, analyzed in Ref. |3lj using normal Gaussian coordinates. The other one is a STT corresponding to the 
conformal coupling which was studied by Noakes 32] using harmonic coordinates. 

Only for the special cases mentioned above the well posedness of the Cauchy problem has been studied in the 
Jordan frame. However, in that frame, none of the analysis performed so far have envisaged the possibility of writing 
the STT in full 3+1 form in order to proof a well formulated Cauchy problem 2 . One exception is perhaps the analysis 
by Sheel, Shapiro & Teukolsky [33, yj]; who studied numerically the BD theory in the Jordan frame by using a 3+1 
decomposition, but this was performed only for the spherically symmetric case with a very particular gauge choice 
and as mentioned for a very specific STT (more comments about this analysis are given in Sec. V). 

In this paper, another alternative is given for treating the Cauchy problem of STT. The STT treated here are of 
general kind and not only a specific case where the non-minimal coupling (NMC) function takes a particular form. 
Moreover, the field equations of STT are written directly in the Jordan (physical frame) and show (contrar y to the 
common believe) that one can obtain a 3+1 formulation (like the Arnowitt-Deser-Misner formulation of GR |35l l36| ) 
leading to first order in time evolution equations for the gravitational and the scalar field (i.e., a well formulated 
Cauchy problem). In fact, this is the first part of a research in which the Cauchy problem of STT is analyzed. In a 
future work, the plan is to obtain a full first-order (both in time and space) formulation of STT in the Jordan frame 
and perform a hyperbolicity analysis using different gauges in order to establish the well-posedness of the Cauchy 
problem (see Appendix A for the well-posedness using a covariant approach). 

In the author's view, the reason for being interested in analyzing the Cauchy problem in the mathematically more 
involved Jordan frame, relies in that such a frame can provide much more physical insight than that the Einstein 



1 Here a distinction is made between a well formulated and a well posed Cauchy problem. Basically a well posed problem not only requires 
that the system of partial differential equations be put as a system of first order in time but also as first order in space and that the 
characteristic matrix associated with the corresponding full first order system satisfy certain requirements as well [see Eq. 1741 and the 
discussion below]. 

2 In the Einstein frame the 3+1 formulation of STT is almost straightforward since one only requires the 3+1 formulation of GR and 
considers the metric there as the non-physical metric. 
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frame (cf. Ref. [34( for a similar point of view). This latter tends sometimes to hide situations where the (true) 
physical metric can have pathologies. Furthermore, when working in the Einstein frame, one can encounter examples 
where one is not able to recover the physical frame due to the presence of singularities |3JJ, |40j or in situations where 
the maps from one frame to the other [cf. Eqs. © and Q)] are ill-defined 3 . Such undesirable situations might 
happen when F(<fi) is not positive definite. Clearly, under a well defined and physically reasonable STT both frames 
should be acceptable and it is a matter of mathematically convenience to use one over the other (cf. Ref. ^^1). In 
this regard, it is to be stressed that it is not the aim of t his pap er to feed the long debate between both frames. This 
has been the issue of countless arguments (see Refs. [2a,|42,|43| for a review). The purpose here is to provide a new 
alternative which gives a new point of departure for researchers interested in numerical analysis of STT or in other 
applications where the 3+1 formulation is better suited than the covariant approach. 

The paper is organized as follows. In Sec. II, scalar-tensor theories of gravity for a single scalar field are introduced 
in several parametrizations (including the Einstein and the Jordan frame representations). Nevertheless, the rest of 
the paper focuses only in the Jordan frame. In Sec. Ill, the standard 3+1 formalism is briefly reviewed and applied 
to the field equations of Sec. II. The new 3+1 equations for STT are then derived. In Sec. IV, several gauge choices 
(lapse and shift) are proposed for a further numerical analysis. In Sec. V, particular examples of astrophysical and 
cosmological interest are discussed. Finally in Sec. VI, several lines of research are suggested as a natural continuation 
of this work. Various appendices complement the paper. Appendix A shows that the Cauchy problem of STT is well 
posed at least in a particular gauge. Appendix B re-analyzes the linear limit of STT using a covariant and a 3+1 
approach. Some aspects of the weak-field limit are considered. Appendix C presents the 3+1 expression of differential 
operators appearing in several gauge choices. Lastly, Appendix D, sketches a self-consistency check performed on the 
3+1 equations of Sec. Ill using the Bianchi identities and therefore the conservation of the effective energy-momentum 
tensor of STT. 



II. SCALAR-TENSOR THEORIES OF GRAVITY 



The general action for STT with a single scalar field is given by 

1 , .. (\ 



S\g ab ,ct>M= j ^Y^/(<l>)R-\^(^^ 2 +V(cP)jjV^(fx + S matt [g ab ,ilj} , (1) 

where "0 represents collectively the matter fields (fields other than </>; units where c = 1 are employed). 

The representation of STT given by Eq. Q is called the Jordan frame representation. One can parametrize the 
same theories as 



where 



$ := F(<j>) , (3) 
W 



<"bd(*) ■= -T^T- > ( 4 ) 



A(*) := -2=^ , (5) 
here ' indicates d$. For instance, the Brans-Dicke theory with w =const. corresponds to F = 2ttGo(/> 2 /oj and V(<fr) = 0. 



3 There exists also the case where the Brans-Dicke parametrization of STT [cf. Eq. together with the use of the Einstein frame can 
lead to ill-defined maps between the Jordan and Einstein frames |41| . 
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It is also customary to parametrize STT in the so-called Einstein frame by introducing non-physical fields as follows, 



Jab 



F{4>)g a b , 



4ttG 



1/2 



W(<j>* 

F*(4>* 

so that the action Eq. JIJ takes the form 

1 



4ttG V* ((/>*) 



F* 



F{4>) 



16nG 



R* - 2(V* 



4W(<t)* 



^Ttfatx + S matt [g* ab /F* (<f ) , V] , 



(G) 
(7) 

(8) 
(9) 

(10) 



where all quantities with '*' are computed with the non-physical metric g* b and 4>* . 

One can remark that although the field equations obtained from the Einstein frame are simpler than those from 
the Jordan frame, in the sense that the field <jf appears to be coupled minimally to the non-physical metric, the 
matter equations derived from the Bianchi identities V*G* ca = will have sources, i.e., V*T^ ca ^ 0, where ab = 
T« b /F* 3 ((j)*) is the non physical energy-momentum tensor of the matter fields i\>. Nevertheless, in the Jordan frame 
the matter equations resulting from the Bianchi identities V c G ca = turn to satisfy V C T™ = 0, reflecting explicitly 
the fulfillment of the Einstein's weak equivalence principle (this is the origin of the name "physical metric"). 

In the following only the Jordan frame representation of STT will be used (for a detailed analysis of STT in the 
Einstein frame see Ref. |l|). The field equations obtained from the action Q are 4 



Gab 



where G ab = R ab - \g a bR and 



Tab 



T 



T 



f 



GcS 
Gq 

V n (/'V b 



8irGoT ab , 

If'R = V , 



rp<p I rrimatt 

2 ab "T" 1 ab 



g ab V c (/'V c 



Tib ■= (V o 0)(V fc 0)-5a6 



(V0) 2 + V{<t>) 



r - J- f- F 

cff ' 8nf 1 1 ■ SttGq 



(11) 
(12) 



(13) 
(14) 
(15) 

(16) 



Using Eq. (|llfl , the Ricci scalar can be expressed in terms of the energy momentum tensor Eq. i|13|) and then Eq. 1121) 
takes the following form, 



D , _ iv - w - \r (i + 3p (v^) 2 + jrz^tt 

where T ma tt stands for the trace of Tj^ tt and the subscript "matt" refers to the matter fields (fields other that <fi). 
Now, the Bianchi identities imply 

V c T ca = . (18) 



4 Latin indices from the first letters of the alphabet a, b, c, ... are four-dimensional and run — 3. Latin indices starting from i (i, j, k, ...) 
are three-dimensional and run 1 — 3. 
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However, the use of the held equations leads, as mentioned before, to the conservation equations for the matter alone 

V c T™ tt = , (19) 

which implies in the case of test particles, that neutral bodies are subject to no other long range forces than the 
gravitational ones (free falling particles). In other words, test particles follow the geodesies of the physical metric and 
the whole effect of the scalar-field is reflected in the modification of the geometry. 

In the next section the field equations will be recasted in a 3+1 form which is specially suited for numerical 
applications. 

III. THE 3+1 FORMULATION 

The aim of this section is to consider the 3+1 or Arnowitt-Deser-Misner (ADM |35j|l formulation of general relativity 
a la York |3fij | and apply a similar formalism to reformulate the scalar tensor theories of gravity in a way as to obtain 
a set of first order time evolution equations and a set of constraint equations for the initial data. 

One then considers a spacetime (M, g ab ) (assumed to be globally hyperbolic) which is to be foliated by a family of 
spacelike hypersurfaces E t parametrized by a global time function t. The foliation is achieved in the following way. 
On a Cauchy surface E t , one is given an initial data set that satisfy some constraint equations, and then the spacetime 
is "reconstructed" by evolving the initial data using a suitable set of evolution equations. This set of constraints and 
evolution equations are known as the ADM equations of GR. In order to find the equivalent set of equations for the 
case of STT a similar algebraic and geometric decomposition of the field equations (|11|) and i|17|l have to be performed. 
In order to do so, some basic aspects of the usual 3+1 decomposition of the Einstein equations are introduced. 

The procedure for obtaining the 3+1 splitting of the Einstein equations consists in projecting tensor fields in the 
direction parallel and orthogonal to the timelike unit vector field n a (n a n a = — 1) which is normal to E t . The 
projection onto E t is performed by defining the projector 

h a b = 5 a b + n a n b . (20) 

The property of this tensor field is that it is idempotent h^h b = h b . A tensor field 3 T aia3 '" ah b ^ b is said to be 
tangent to £( if contracted with n a is zero or if contracted with h a b remains unchanged. For brevity, this tensors will 
be termed as 3— tensors. Any tensor field can be decomposed orthogonally using h\ and n a . In particular, a 4- vector 
w a € T^ 1 is decomposed as follows 

w a = 3 w a + w ^ n a ^ ^ 

where 3 w a :— h a b w b and w± := —n c w c , with signg ah = (—,+,+, +). Moreover, the 3+1 splitting of the metric is as 
follows: 

ds 2 — — {N 2 — AT*iVi)d£ a - 2N l dtdx l + h lj dx l dx j , (22) 

where the lapse function N > is defined as to normalize the (future pointing) dual vector field n a = —NV a t- The 
shift vector is given by N a := —h a b t b , where t a is a vector field that represents the "flow" of the time lines and that 
satisfies t a V a t = 1. This means that t a is orthogonally decomposed as t a — —N a + Nn a , with N — ~n a t a . Here hij 
is the 3- metric (or induced metric) of the manifold S t . 

Another important object is the extrinsic curvature of the embedings Et which is defined as 

K a b '•= —-^C-nhab , (23) 
where Cn stands for the Lie derivative along n a . From the above definition one can obtain the following identity 

K ab = -h a c h b d V c n d , (24) 

which shows that K ab is a 3-tensor field. 

As it is well known, the set (E t , h ab , K ab ) provides the initial data for the gravitational field. This data cannot be 
arbitrary but have to satisfy the Einstein constraint equations (see below). 
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At this point it is useful to introduce a derivative operator compatible with h a \,. Given a 3-tensor field 
3 7 iaia2 - a W.. 6i , one defines 

«e 3 ^-" o Vi a ...6, = ^^»^C^-Af , ^v / 8 r ^-" c \^ , (25) 

where D c h c d = 0. Finally, it is to be mentioned that the indices of 3-tensors can be raised and lowered with h ab and 

h a b- 

For the purpose of this paper, it is convenient to introduce the following quantities 

Q a := D a 4> , (26) 
n := C n <t> , (27) 

where IT is to be identified as the momentum of the scalar field which plays a similar role as K a b for the gravitational 
field. On the other hand Q a is the analogous of the 3-Christoffel symbols in the sense that it contains first order 
spatial derivatives of the scalar-field. Nonetheless, unlike the 3-Christoffel symbols, Q a is as a 3- vector field. 
Now from Eqs. I|23[l and (|27|l one obtains the following useful expressions 

KlJ = -V inj = -JVrfc = + D 3 N % + AiV,) , (28) 

n = n a V a 4> = 1 (dtt + N-Qa^j . (29) 

Eq. (|28|l is to be regarded as an evolution equation for while l|29|) provides the evolution equation for <\>. In fact, 
one can obtain an evolution equation for Qi by applying Di on l|29|l and using df t 4> = d^<p = dtQi, where in this case, 
Di = h1d c = di. Therefore 

d t Qi = Di(NIL) - Di(N l Qi) . (30) 
Because of the symmetry DiQj = DjQi which follows from the integrability condition dfjtfi = dj^, one can write 

d t Qi + N l diQi + QidiN 1 = Di(NU) . (31) 

This is indeed equivalent to 

CnQa = ^D a (NU) . (32) 
The orthogonal decomposition of the energy- momentum tensor is given by |3fij | : 

T ab = S ab + J a n b + n a J b + En a n b , (33) 



where 



S ab := h a c h b d T cd , (34) 



J" 
E 



c d 

' a rricd n 



= -h\T cd n d , (35) 
= T cd n c n d , (36) 



here S ab is the 3-energy-momentum tensor, J a is the ^-momentum density vector and E is the total energy density 
measured by the observer orthogonal to S t with 4- velocity n a . 

In this paper T ab is the effective energy-momentum tensor (II 3D which includes the contribution of the scalar field 
and another kind of matter: T ab = J2i ^7 b - This means that 



i i i 



(37) 



7 



Specifically for the analysis at hand the above quantities are given by 



E 

Ja 
Sab 



^(Ef 



Go 
Go 

G c ff / / 



Jf- 



S 



ab 



E r - 



' it 



oniatt \ 
dab J 



(38) 
(39) 
(40) 



In order to proceed with the 3+1 formulation of STT, one is required to project orthogonally the Einstein equations 
ffTTf) - The projection of Einstein equations G a b = 8nGoT a b [or R a b — 47rGo (2T a b — T c c g a b), better suited to obtain 
the dynamic equations] in the directions tangent and orthogonal to S t , followed by the use of the Gauss-Codazzi 
equations leads to the 3+1 form of Einstein equations a la York |36[ : 



known as the Hamiltonian constraint. 



known as the momentum constraints. 
The dynamic Einstein equations read 5 



d t K\+N l diK l 



R + K 



DiK 1 



K\d N l 



l6irG Q E 



DiK = %kG J 1 



(41) 



(42) 



K l ,,diN l 



3 rti 



R l jN — NKK l j 



= 4nG N 



(S - E)8 l 1 - 2S\ 



(45) 



where K and S stand for the traces of the extrinsic curvature and the 3-energy-momentum tensor respectively. All 
the quantities written with a '3' index refer to those computed with the three- metric hij. 

An evolution equation for the trace K is obtained by taking the trace in Eq. (|43|) followed by the use of Eq. l|4T|) 
to give 



d t K + N l diK 



AN - NKijK 13 



AttGoN 



S + E 



(46) 



where 3 A := D l Di stands for the Laplacian operator compatible with h 



As mentioned at the beginning of this section, in pure and vacuum GR, given an initial data set (T, t ,hij, Kij) 
satisfying the constraints Eqs. I|41|) and 142|) . one performs the Cauchy development by evolving hij and K{j forward 
in time using Eqs. l|28|l and l|45(l . respectively. Of course, in order to achieve the evolution, lapse and shift conditions 
(gauge conditions) have to be imposed (see Sec. IV). 

Now, although the Einstein equations give in general no information about the specific form of the energy-momentum 
tensor (EMT), one usually considers some features on T a j, that leads to a well posed initial value problem. In particular 
is common to consider EMT that does not contain second order time derivatives (SOTD) of the matter fields so that 



5 This can be written as 



or alternatively as 



h a c h b d C n K c d + jfD a D b N - 3 R a b - KK% = 4ttG 



(5 - E)h\ - 2S a b 



C n K ab + -D„ D b N - 3 R ab - KK ab + 2K ac K\ = AnG 



(S — E )h ab - 2S ab 

where C n K ab = {C t K ab + C N K ab ) or h a c h b d C n K a d = jfh a c h b d {C t K^ d + C N K c d ) 



(43) 
(44) 
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the theorems for existence and uniqueness of solutions of Einstein equations together with the matter equations 
V a T ab = can be applied (e.g. see Refs. [HHU ). 

Now, in the case of STT, the non-minimal coupling (NMC) of the scalar-field with the curvature makes however 
that the right-hand-side (r.h.s) of Eq. actually contain second order derivatives of (f> both in time and space. 

Therefore, a well formulated initial value problem seems in jeopardy. In order to obtain a nice initial value formulation, 
such terms (notably, the SOTD terms) should be eliminated using the remaining field equations. In what follows 
such elimination is performed [namely, in Eq. (jl4|l ] in order to obtain the explicit form of the 3+1-effective-EMT. 
The remaining contributions of the EMT Eq. 1|13|) are quite standard and so their 3+1 decomposition pose no big 
problem. In this way the goal is to isolate the pure matter terms and the scalar-field sources at the r.h.s of Eqs. 
(|41|l — (|45|l . and then to place the other quantities associated with the scalar field in the left-hand-side (l.h.s.) so that 
the scalar field be treated more less at the same footing as the gravitational field. 

So let us proceed to compute E$ — n a n b T^ b . A straightforward but lengthy calculation leads to 



E* = /' ( 3 A0 + Kn a V a <t> ) + f"(D(j)) 



where {D<j)) 2 := (D c <t>){D c (j)). Note that E f no longer contains SOTD. 

Now one computes Sj b = h a c h b d Tj d . An equally lengthy computation yields 



oab 



f [l3 a D b cj) + K ab (n c \7 c( j)) - h ab Ucf^ - f" {h ab [(L>0) 2 - (n c V c 0) 2 ] - (£>»(£>»} 



(47) 



(48) 



As one can appreciate, SOTD appear in Sp b only through □ </>. However, by using the r.h.s of Eq. I|17|) one can replace 
the d'Alambertian in terms of first-order derivatives and sources. This is done below. On the other hand, the term 
D a D b (j) contains only second order spatial derivatives (SOSD). 
From the above expression one can easily calculate its trace 



S f = f 



! A0 + K{n c S7 c cj>) - 3D0 



3(n a V a ^) 2 - 2(D<f) 2 

1 (n a v„# 



From H47|) and (|49|l one can obtain the quantity 

S f - E* = -3/'D0 - 3/" \{D<t>) 
Finally, one requires J/ = —n d h^T^ d . This quantity reads 

J[ = -f [K a c D c <t> + D a {n c V c 4>)] - /"(n c V c 0) D a 
Using Eqs. and in Eqs. I07J)-(|3TJ one obtains respectively 
E f = f [D a Q a + KTl) + f'Q 2 , 

/' (D a Q b + IlK ab - h ab Ocj)) - f" \h ab (Q 



Sf 



n 2 



Q a Q l 



S f = f [D a Q a + KU- 3D0J + /" [3n 2 - 2Q 2 

E' = -3/'n^-3/"(Q 2 -n 2 ) , 



where Q 2 := Q C Q C . 

In a similar and straightforward way one finds 

E,h 



nab 

D 4, 



2 

Q a Q b 



n 2 + Q'- 



V{cj>) , 



h ab 



n 2 ) + v(4>) 



S,p 

Scf, — Eff, 



\ [m 2 - q 2 ) - 3V(4>) 
n 2 - q 2 - 4V(4>) , 

-TlQa ■ 



(49) 
(50) 

(51) 

(52) 
(53) 
(54) 
(55) 
(56) 

(57) 
(58) 

(59) 

(60) 
(61) 
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Now, collecting these results, and according to Eqs. I|38f) — (|40|l and using Eqs. I|l(jfl and l|17|) . the total contributions 
of the 3+1 effective energy- momentum tensor read 



E 

Jo 



1 



8ttG / 
1 

8ttG / 
1 



/' [D C Q C + + + 2/") + V{<£) + E n 



- f (k;q c + D a n) - uQ a (i + /") + j a matl 



(62) 
(63) 



Sab = { QaQb (l + /") + /' [DaQb + IlK ab ) - 



h a b 



3/1 2 
2/ 



-n 9 ) (1 + ^ + 2/ 



+v ( 1 _ i? ) + /V + 4? ( 5matt _ ^ matt 



2/ 



2/ 



crnatt 
D ab 



It turns useful to write explicitly the following quantities because they appear in different applications 



S-E = 



8ttGo/ (l + ?ff- 

s + E = ^ f { 2f i D c QC + UK ) + J^E 



n 2 - Q 2 ) (l + 3/") - 3fV - W{4>) + 5 matt - E matt 



(64) 



(65) 



2IT 1 



*f'\ V" U 9H2 ( 
4/ + 2 +2C M 4/ 



f" 

(1 + 2/")- V 



3f 



3r 



-2V I 1 - -±y- J - 3f'V + SUtt + ^matt [ 1 + ~J 



(66) 



The input of these expressions in Eqs. (|41|) . (|42|l . H45|) yields finally the 3+1 equations of the scalar-tensor theories 
of gravity given by the action (Q. The Hamiltonian constraint, the momentum constraints and the dynamic Einstein 
lr T are resnectivelv. 



equations of STT are respectively, 

3 R + K 2 - K tl K l] - - 
f 



f ( A Q i + /m) + ^ + ^(i + 2/ 

/' (KjQt+DJl) +UQ i (l + f" 



DiK 1 i — DiK + — 



K 



matt 



-v{4>) 



d t K' l + N l diK'j + K\d 3 N l - K l &iN l + D i D j N ~ 3 R* 3 N - NKK 1 



N 

T 



Q l Qj (i + f")+f (d'Qj+uk 1 



3/ /2 
2/ 



2/(l 



3£! 
2/ 



1 rmatt 

"7 1 



J ill 



(67) 
(68) 



N 



25 



matt j 



f'V +2V il+ J — \ - (S matt - £mat, 



r 
/ 



.(69) 



On the other hand, Eq. I|46|l reads 



d t K + N l diK + 3 AN - NK i3 K lj - ^f- (DiQ 1 + UK 



f 



N 



'n 2 1 



3/' 2 , 3/" 



S 



'matt H~ ^matt 



3T 
4/ 



(l + 2/" 



/ 



■2V 1 



3T 
2/ 



r 
2 



3/V 



(70) 
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Finally, Eq. (|17fl takes the form (c.f. Appendix C) 

fV' - 2fV - if (1 + 3/") (q 2 - n 2 ) + i/'T K 
CnH -IIK- Q c D c [lnN] - D C Q C = -. '- , (71) 



2/ 

where £ n II = -^<9 f lT + ^<9;1T and D C Q C = d c Q c + Q c d c (lnVh), here h = det/ij, and T matt = 5 ma tt - E matt . 

The Eqs. (|29|) . (|31|1 and l|71|l are the evolution equations for the scalar-field variables. In fact, since one has 
introduced an evolution Eq. (|32|1 for Q a , formally one can also promote the definition of Q a Eq. Ij26(l and the 
integrability condition D a Qf, = DbQ a as constraints. Taking the spatial components one has 

Qi-Di<t) = , (72) 
D[iQj] = . (73) 

The system of evolution equations (both for the gravitational field and the scalar field) are to be completed with 
the evolution equations for the matter. 



In summary, the Cauchy problem of STT can be sketched as follows: given the initial data set (E t , hij, Kij, <f>, II, Qi) 
subject to the constraints Eqs. lE3, (JBSJ, and tfTHft . the set of equations P5). (JUJ), lEH), EH and JHJ provide 
the Cauchy development for the gravitational and the scalar field. Clearly this development can be unambiguously 
performed once the gauge (choice of coordinates) has been fixed. This is equivalent to fix the lapse and the shift. 
Note that these quantities appear in all the evolution equations. In the next section, the gauge choice for STT is 
proposed by generalizing several lapse and shift evolution equations that have been analyzed in the past. 

Now, more comments about this new set of 3+1 equations for STT are in order. A nice feature of the constraint 
Eqs. (ff)?j) and ljf)%)l [see also fT2|l and (f73^l ] is that the lapse and the shift do not appear explicitly. These equations 
then constraint the initial data hij, Kij, (f>, II, Qi) regardless of the coordinate choice. On the other hand, neither 
the constraints l|67|) and (|68|l [see also l|72|l and 1|73[) ] nor the dynamic Eq. I|69|) contain time derivatives of II or 
Q a - Moreover, the system of 3+1 equations [notably the dynamic Eqs. i|t)9[) and l|71l) ] are now of full first order 
(in time and space) in the scalar-field variables. This is to be contrasted with the pure gravitational sector of the 
equations where SOSD of hij appear when expanding the curvature 3 R l j or in the explicit term D l DjN for the lapse. 
More recently, new 3+1 formulations (e.g. see |44ll45| and references therein), provide a full first order system of the 
gravitational and gauge sectors which are better suited for numerical stability (see the discussion below). 

In order to avoid any kind of pathologies one should only consider STT with / > 0, particularly if one requires 
a mathematically well posed initial value problem (see Appendix A). For instance, the condition / > avoids the 
possibility for the appearance of potential singularities in Eqs. I|67|) . H68[) and l|tj9|) . Notably in the terms with 1// and 

1/ ^1 + 3f' 2 /(2f)j. Those kind of singularities could preclude the existence of a Cauchy hypersurface 6 . Moreover, 

physically the condition / > ensures a positive definite effective gravitational constant G e g [cf. Eq. (|l(j|) ]. 

It is to note that if / = const, the new set of 3+1 equations reduce to the standard 3+1 equations of GR with 
a minimally coupled scalar-field. On the other hand, as discussed in Appendix D, the full form of 3+1 equations of 
STT have passed successfully a self-consistency check that consists in computing the 3+1 conservation equations of 
the effective energy-momentum tensor of STT and showing that such equations together with the 3+1 field equations 
lead to the conservation of the energy-momentum tensor of matter alone. 



For instance, an initial data set corresponding to a static and spherically symmetric spacetime with / = (1 — 47rGo<ji 2 /3)/(87rGo) and 
V(<fi) = 0, and without matter, is the one associated with the well known solution of Bekenstein-Bronnikov-Melnikov-Bocharova |4(I . 
where the spacetime has the geometry of the extreme Reissner-Nordstrom solution. In this solution, the scalar field is singular at the 
event h orizon. Moreover, it has been argued that the Einstein field equations are not satisfied at the event horizon l60l (see also Refs. 
39, 4C| and references therein for more considerations about this point and the conformal coupling). Those kind of pathologies could 
be avoided with the inc lusion of matter (e.g. a cosmological constant A with V((j>) = A</> 4 and A = — 27rA/9 can give rise to regular 
black- hole solutions 1471 ). 
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The set of evolution equations (|28|) and 145|) in the absence of matter fields and the scalar field (i.e, vacuum GR) 
cannot be straightforwardly submitted to an hyperbolicity analysis which is of main importance for establishing a 
well posed Cauchy problem. Only recently the 3+1 equations (including evolution equations for the gauge functions 
-see Sec. IV-) have been rewritten in full first order form, 

d t u + U i V i u = S(u) , (74) 

where u represents collectively the fundamental variables (like the hij's, Kij's, etc.), MP is called the characteristic 
matrix of the system and S(u) represents source terms which include only the fundamental variables (not their 
derivatives). This system of partial differential equations (PDE) is said to be quasi-linear since it is linear in the 
derivatives but in general it is non-linear in u. A system of PDE of the form Eq. I|74ll is said to be: 1) weakly 
hyperbolic, 2) strongly hyperbolic 3) symmetric hyperbolic or 4) symmetrizable hyperbolic if the characteristic matrix 
MP corresponding respectively to such cases satisfies: 1) Given any co- vector Si, then SiM 1 has a real set of eigenvalues 
but not a complete set of eigenvectors, 2) SjM 1 has a real set of eigenvalues and a complete set of eigenvectors (plus 
some bound conditions), 3) IMP are symmetric, 4) M 1 can be symmetrized. 

Only strongly, symmetric and symmetrizable hyperbolic formulations (hereafter hyperbolic formulations) have a 
well posed Cauchy problem (see Ref. |28| for a review). On the other hand, quasi- linear wave equations (QLWE) 
are hyperbolic and they have a well posed Cauchy problem |49|]. The so called second order analysis (as opposed to 
the first order one outlined above) consists in re-writing the evolution equations as QLWE. For instance, as it is well 
known, the vacuum Einstein's field equations reduce to a system of QLWE in the harmonic gauge. 

After many years of using the 3+1-ADM equations as a "work-horse" for numerical studies, now it is known that 
in general Eqs. (|28[1 and l|45|l (in vacuum) even in its first order form like in Eq. (|74|) . are not hyperbolic (only weakly 
hyperbolic) and that one needs to add some multiples of the constraints to turn the whole system hyperbolic or 
at least to obtain a numerically stable system (see Ref. ^^). Of course this also requires the use of suitable gauge 
choices (see Sec. IV). In this sense, it is expected that the new system of evolutions equations (|28ll . I|29(l . (|31|l . 1691) 
and (|71|l [once put in first order form J73J] will not be a hyperbolic system either. Clearly the pure scalar-field sector 
is hyperbolic since Eq. I|17f) is a quasi-linear wave equation. 

It is out of the scope of the present work to provide a full first order formulation of Eqs. (|28ll . (|29l) . 131|) . I|69l) and 
(|71|l and its corresponding hyperbolicity analysis (e.g., see Refs. for the case of GR). Nonetheless, the new set 

of 3+1 equations is a first step towards this task. On the other hand, using the the covariant equations l|ll|l and i|17fl . 
together with a generalization of the harmonic gauge for STT, it is shown in Appendix A that STT in the absence of 
matter have a well posed Cauchy problem. 



IV. SLICING AND SHIFT CONDITIONS 



As it was mentioned in previous section, the evolution of the gravitational and scalar fields can be performed 
without ambiguity once the lapse and the shift are fixed. This fixing amounts to impose coordinate conditions. 

Let us remind that one of the most convenient coordinate (gauge) choices in GR which leads to a well posed initial 
value problem is the De Bonder or harmonic gauge (HG) given by 

□/ = o = -La c (V^r) = -rt . (75) 

In the framework of the 3+1 formalism where the coordinates x a are chosen to coincide with (t,x l ), that is, with 
the time global function that defines the spatial slices S t and the local spatial coordinates x l of the submanifold E t , 
the HG leads to an evolution equation for the lapse N and the shift iV a (see below). For the analysis of the Cauchy 
problem of STT, a generalization of the harmonic gauge is proposed [hereafter pseudo-harmonic gauge (PHG)] in the 
following way 

Ux a = ~Q\7 a T , (76) 

where T := InF and the r.h.s of Eq. (|76|l is to be understood in the following way V a JF = g ab Wb^F- Here O is 
an arbitrary parameter. As analyzed in the Appendix A, for = 1 this gauge leads in fact to a well posed initial 
value problem of STT. Moreover, in the linear-limit (see Appendix B) this gauge choice (with = 1) generalizes the 
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Lorentz-De Donder like gauge which loads to wave equations for the spin-2 and spin-0 modes (gravitational and scalar 
waves propagating in "Minkowski" spacetime). 

In numerical relativity, it is customary to chose the lapse and the shift conditions independently. This means that 
one does not necessarily impose the same condition Eq. (|75|l simultaneously for the time and spatial coordinates. It 
is in fact sometimes better to impose rather different equations for both the lapse (slicing condition) and the shift. 
One of the most popular family of slicings is the so-called Bona-Masso (BM) conditions [5l| which can be written in 
four-dimensional notation as [52j : 

at = T q n b n c VbV c t , (77) 

where T q = 1/q—l and q = q(N) > (i.e., q is a positive but otherwise arbitrary function of the lapse). This BM 
condition provides the following evolution equation for the lapse (see Appendix C) 

1 N l 

CnN = -d t N + —diN = —Nq(N)K . (78) 

More recently, it was proposed a modification of the BM condition in the following way 13] : 

T /N c \ 
□i^tVViVci-Ve^J . (79) 

This modified Bona-Masso (MBM) condition is equivalent to 



d t N = —N q(N) \ K + —D c N c j . (80) 

This condition has the advantage that the lapse does not evolve if there exists a timelike Killing vector field (d/dt) a . 
In that case dthij = 0, [in particular for the determinant dth = —2h (NK + D C N C ) = -cf. Eq. H28fl - ] and therefore 
the r.h.s of Eq. 180|) vanishes identically. 

For the BM or MBM, different choices of the function q(N) have been analyzed both analytically and numerically 

PE10. 

For the case of STT, the following lapse and shift conditions are considered (from which several gauges in GR can 
be recovered): 



Ux a = Ar qa n b n c V b V c x a +BS% 



T / N c \ N c 



-8V°J (no summation over the index a) 



. (81) 

Here one allows to have four different functions T qa — l/q a — 1, with q a = q a (N) > 0, and where A and B are arbitrary 
coefficients. For instance, if B — and A — 1, one has a STT generalization of the BM condition for x° = t. If A = 
and B = 1 one finds the STT generalization of the MBM condition for x° = t. The pure GR gauge conditions are 
recovered trivially when T — const. 

A straightforward but cumbersome algebraic manipulations using the 3+1 formalism of Sec. Ill (see Appendix C), 
allows one to re-write Eq. H81|) as the following evolution equation: 



T ^ { dcna + K ") + ^ ( DcNC - NCD c(^N) 



nW-D>Q*N) = 1 + T J A + m) {v a + n«K + B5 ? 

+Q^(Q--n a U + S^N c Q ( }j^ , (82) 

where one reminds n a — (1/N, N l /N). Here V a :— -^d c (^Vhh ac ^j, and it was used the following identity V a JF = 

T'\J a cf> = T' (D a <p - n a n c V c <p) = T' (Q a - n a Il). Note for the time "component" V* = 0. 

The time and the spatial parts of Eq. Ij82(l yield respectively the evolution equations for the lapse and shift: 

W + /.V^Tm ^gj* = - 1 I ^ , m W + ^ f n - \ > (83) 



l + T qt (A + B) J l + T qt (A + B) { N L v N 
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n c d r I — ] - D l {\nN) = - - \ V 1 + —K + QT' [ Q l - —U 



N 



1 + T gj ^l 



N 



N 



(84) 



When the lapse and the shift equations (|83|l and (|84|l are imposed simultaneously and for the simple case where B = 
and the four functions T qa = T q are all equal then Eq. I|84|) can be simplified using Eq. I|83|l to yield 



d t N l + Wd.N 1 - ND i N = 



N 



\v 1 + qf'q 1 } 



which is valid if the lapse satisfies 



d t N + N { diN 



N 2 



{k - e^'n} 



It is also usual to employ the variable N % = N l /N instead of the shift {55|. Then, Eq. (|84|l reads 

N 



d t N l + NN ] djN l - D l N 



1 + T 9i A 



V 1 + N l K + Of Q l - N'U 



(85) 



(86) 



(87) 



Again, note from Eq. I|86|l if A = 1 and T q = l/q(N) — 1, one has a generalization of the BM slicing family for 
STT: 



d t N + N l diN = —q(N)N 2 (^K - QT'TY 



In the same way, note from Eqs. H85JI an d l|87|l that if A = 1 and T q = l/q(N) — 1 and T qi — l/p(N) — 1 one 
obtains the generalization for STT of the shift evolution equations proposed in Ref. [55| : 



d t N l + N j djN l - ND i N = q(N)N 2 (v* + QT'Q^ , 



(89) 



d t N l + NN j d 3 N l - D l N = p(N)N 



V 1 + N l K 



-N l n 



(90) 



On the other hand, note also from Eq. I|83l) that for the case A = and B = 1 and T qt — l/q(N) — 1, one obtains 
a generalization of the MBM condition for STT: 



d t N = ~q{N)N 2 



1 / N c 

K+-D C N C -QT' \ - —Q, 



(91) 



which, in the presence of a timelike Killing vector field (d/dt) a , also leads to a time-independent lapse (like in the 
original MBM condition) because of the arguments given before and also because in that case d t (j)=0 which by virtue 
of Eq. HU) then leads to TVn _ atcq c = fj. 

Finally, if q(N) — 1, then from Eqs. (|88|) and 189|) one recovers the PHG Eq. Q76f) here in terms of the lapse and 
shift 



d t N + N l diN = -N 2 - QT'Uj , 

d t N l + N j djN' - ND i N = N 2 (v l + QT'Q^ . 



(92) 
(93) 



In Appendix B the linear limit of Eqs. (|92|l and (|93|l will be taken and show that this pair of equations generalizes 
the Lorentz-De Donder gauge of GR to the case of STT. In pure GR {J- = const.), the set of evolution Eqs. I|92|) and 
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(|93[l reduce to the harmonic gauge that was considered in the past by Smarr & York [5(j in terms of the lapse and 
shift 7 . 

In the case of GR, the BM and MBM conditions (given a prescribed shift) have been shown to be consistent with 
hyperbolic formulations [5ll Isif . Alternatives to the slicing conditions are the use of a densitized lapse S = ln^/i"' 3 ] 
|44| (where 2(t,x l ) is an arbitrary but a priori known function of the spacetime coordinates and (3 is a constant; a 
positive f3 is favored for strong hyperbolicity) or to take the lapse N(h, t, x l ) as a prescribed function of the determinant 
h of the 3- metric and the spacetime coordinates 56] (dhN > is required for strong hyperbolicity). 

On the other hand, it is important to emphasize that in pure GR (not even to mention the STT) only a few shift 
evolution equations have been tested numerically in the past or included in the analysis of hyperbolicity |57l l58l l59| . 
A particularly simple shift is a null one (i.e., N l — 0, d^N 1 = 0; gauge not included as a special case of the above 
shift conditions unless one considers a constrained initial lapse). However, this shift choice is not always free of 
numerical pathologies. Shift conditions have not been explored thoroughly, and only empirical choices have proved to 
be successful in numerical experiments. The shift equation <|9(Jll in GR (J- = const.) have been analyzed numerically 
for very special cases [55|. 

There also exists gauges leading to elliptic equations rather than wave equations for N and N l . One is the maximal 
slicing condition (K = 0, dtK = 0) which in GR would correspond to take q(N) arbitrarily large [cf. Eq. I|78ll]. Then 
Eq. H46fl provides an elliptic equation for N. An elliptic equation for N l was considered by Smarr & York [50J (which 
includes the so-called minimal distortion gauge) . Both elliptic equations were shown to be the analogous in GR of the 
radiation gauge in electrodynamics. It would be interesting to generalize those conditions in STT. For, in Eq. H9()(l 
taking q(N) arbitrarily large would lead to K = QT'Ii which for 8 ^ and in view of the evolution Eqs. i|7U|) and 
H7H one would obtain a complicated elliptic equation for N. 

In any case, a full and new numerical exploration of the above gauges (or others) is to be performed in order to 
test their usefulness in the context of STT. 



V. SPECIFIC APPLICATIONS 



The new system of 3+1 equations presented in Sec. Ill when applied to specific spacetimes coincides with the 
system of equations analyzed in the past 0, ITEl flrl El) and which were derived from the covariant field equations in 
a rather circuitous way in order to recast them into first order form. Below, some details about some of these specific 
situations are discussed. 



Cosmology. In the case of an isotropic and homogeneous spacetime (Friedmann-Robertson- Walker cosmology) 
with a combination of two non- interacting perfect fluids (representing photons and cold matter), the 3+1 equations of 
Sec. Ill reduce to those of Refs. [HEl for the class of STT where F(<f>) = 1 + 16ttGo^ 2 [£ > 0; recall /(</>) = F/(8tt)] 
and V(4>) = m 2 (j) 2 . In particular, the Hamiltonian constraint Eq. I|t)7|) leads to the Friedmann equation relating the 
Hubble parameter with the total energy-density of the Universe. For this case, the extrinsic curvature is related to the 
Hubble expansion as K = —3a/a [a{t) being the scale-factor] and 3 R — 6k/a 2 (where k = ±1,0). On the other hand, 
the dynamic equations (|69[) lead to a single equation for the acceleration of the Universe (a) . As shown in Refs. [ITl[l6| . 
the two first order equations [one for a resulting from Eq. (|28[) ] and the other for K [resulting from Eq. (7UJ)] together 
with the Eqs. (|29|l and (|71|l for the scalar- field can be integrated numerically. The momentum constraint Eq. (JBSJ, as 
well as the constraints Eqs. I|72l) and (|73|l . are satisfied trivially in this case. The Hamiltonian constraint can be used 
at every time step to monitor the numerical accuracy of the results 8 . Regarding the matter, the matter conservation 
equations together with the equation of state (EOS) for photons and dust lead to the standard quadrature for the 
total energy-density of matter /? ma tt = C\ /a + ci /a 3 (where c\ , ci are constants indicating the content of photons 
and "dust" -baryons and dark matter- at a given cosmic time). This set of equations were successful for analyzing 
the evolution of the Universe by taking as "initial" conditions the conditions of the Universe today and integrating 



7 One should notice a sign error in the second of Eqs. (6.4) for the shift in the Smarr & York article |50H (it should read "— " instead of 
'+" in front of the term 7 ! -T* fc ). For the purpose of comparison with that article or others which use a different sign convention and 

notation for the shift, one should take into account the following N l = —/3 t . 

8 Here it is important to mention a typographical mistake committed in Eqs. (34) and (A19) of Ref. M. Salgado, D. Sudarsky, and H. 
Qucvedo, Phys. Rev. D 53, 6771 (1996) that concerns the sign of the Universe's spatial curvature given in terms of the total-energy 
density. However, the correct sign was taken into account in the numerical code and hence the results are not affected. 
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the equations backwards in time. The initial conditions for the scalar-field and the value of the NMC constant £ can 
be chosen so as to satisfy several local constraints |l6|] . The concrete STT that results is completely rigid (falsifyable) 
in the sense that all the parameters are fixed for once and for all and that the specific STT should account for all 
the observations (which in fact give raise to a larger number of algebraic constraints than the available parameters of 
the theory) . Among these observations one can mention the primordial nucleosynthesis (specially that of He 4 which 
restricts severely the expansion rate of the Universe at early times), the type-I supernova distance-relation, the Cosmic 
Microwave Background measurements, etc. 

At this point it is worth mentioning that a remarkable number of investigations in cosmology using the Brans-Dicke 
theory or other STT exist in the literature. It is completely out of the scope of this work to comment about each of 
them. Here the aim is simply to make contact between the equations of Sec. Ill and their use in cosmology. 

Static and spherically symmetric spacetimes. In the case of a spacetime admitting a static time-like Killing 
vector field and 5*0(3) symmetry, and also for the class of STT with F(<f>) = 1 + WnGo^ip 2 (£ > 0), the 3+1 equations 
of Sec. Ill reduce to those of Refs. [lUl41j . In those references, for the matter it was taken a perfect-fluid representing 
a neutron star. In this case, one is basically dealing with solving the constraints (the Hamiltonian constraint), the 
equation for the lapse (the norm of the time-like Killing vector field) and the Eqs. (|26|l and l|71|) for the scalar-field. 
In this case the quantity DiQ 1 which appears in the Hamiltonian constraint Eq. i|67[l can be replaced in terms of 
simpler quantities using the static limit of Eq. I|71|l . The boundary conditions imposed in Refs. [U|41| were such as 
to obtain an asymptotically flat spacetime and fields regular at the origin. This analysis led to the recovery of the 
phenomenon of spontaneous scalarization discovered by Damour & Esposito-Farese 0, but now using realistic EOS 
and in the Jordan frame. Moreover, a further analysis along this line can be found in Ref. |4lj| (this deals mainly 
with the energy-conditions). 

Within the context of static and spherically symmetric spacetimes, it is perhaps worth mentioning some of the 
situation of black holes in STT. In the early seventies it was somehow a surprise to find that for the conformal coupling 
F(<fi) = 1 — 47rGo0 2 /3 and V((j>) = (and without matter) the field equations admitted an exact asymptotically 
flat static and spherically symmetric (AFSSS) solution (the Bekenstein-Bronnikov-Melnikov-Bocharova solution) jig 
which corresponds actually to the extreme Reissner-Nordstrom geometry, except that the scalar-field replaces the 
electric field. Despite the pathologies of the scalar field at the horizon (it is singular there), it was considered as 
a genuine solution, until very recently. Sudarsky & Zannias |60| carefully reanalyzed the issue and showed that in 
fact due to this pathology the "solution" does not hold at the horizon. The lack of regularity of the scalar field at 
the horizon (a standard requirement for the existence and uniqueness of solutions) was then the responsible of this 
anomaly (see also Ref. 61]). Actually, the regularity at the horizon is also a crucial condition for establishing the 
opposite situation: the non-existence of scalar- hair (the so-called no hair theorems) These theorems apply for 

STT with a minimally coupled scalar-field and non -neg ative potentials V(4>) > for AFSSS. For the case of NMC, 
no-hair theorems also exist for certain class of STT |63| • 

When the condition V(4>) > is abandoned, hairy back holes have shown (at least numerically) to exist |64j . 
Moreover, by introducing a cosmological constant, it is possible to find scalar- hairy black holes in asymptotically 
de-Sitter and anti-de Sitter spacetimes with scalar-fields coupled both minimally and no-minimally |47l l65|. 

Nevertheless, for asymptotically flat spacetimes and with V((j>) > or V(4>) = 0, there is still many open questions 
about the existence of regular black hole solutions in STT. 

In many aspects the new 3+1 equations of Sec. Ill can be useful for the analysis of stationary and static spacetimes. 
In those cases, many degrees of freedom cancel out and the equations simplify considerably. For instance, in the case 
of static spacetimes, ify = = n. The momentum constraint Eq. (|68[) reduces to = 0. In static and spherical 
symmetric spacetimes the simplification is even greater since one can chose N l = 0. In fact, one has only to solve an 
elliptic equation for the lapse [cf. Eq. I|7UI) ; or alternatively a first-order equation for the lapse which results from the 
term D l DjN when using the angular components of the dynamic Eq. I|69l) ]. the two first-order equations for the scalar 
field [cf. Eqs. I|26|) and J7TJ], and the Hamiltonian constraint for one metric potential (e.g., g rr )- Moreover, when 
employed the parametrization g rr — (1 — 2m(r)/r)~ 1 , the Hamiltonian constraint provides a very simple ordinary 
differential equation for m(r). Indeed and as commented above, when taking into account a perfect fluid, it is such a 
set of equations that were used in Refs. [ITll4l| to treat the phenomenon of spontaneous scalarization. 

Dynamical spacetimes. The dynamics of spacetimes in STT have been analyzed in the past only under 
very symmetric cases. In Brans-Dicke theory the spherically symmetric collapse of dust to a black-hole formation 
(Oppenhcimer-Snyder collapse) has been analyzed by Sheel, Shapiro & Teukolsky |33,|3J|. These are as far as the 
author is aware, one of the few numerical studies where a STT is analyzed in the Jordan frame. A 3+1 decomposition 
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of the equations are considered (but only for the special case of spherical symmetry) . The slicing condition employed 
there was maximal slicing, and the spatial coordinates are of isotropic type (with a no-null shift). 

A similar study performed a little bit earlier was carried out by Shibata, Nakao and Nakamura 3 (see also p|), 
also in BD theory and for the Oppenheimer-Snyder collapse. This study, however, used the Einstein frame for which 
the standard 3+1 formalism was applied. The maximal slicing condition was adopted together with a null shift. The 
scalar-wave signal resulting from the dynamics was confronted with the sensibilities of interferometric detectors like 
LIGO. 

More recently, Novak |Tcj performed several studies in spherically symmetric neutron-star models within certain 
class of STT. In Novak's analysis the Einstein frame was used and the corresponding 3+1 decomposition (which is the 
standard one) was implemented for this particular case. The coordinates used there correspond to the simple radial- 
polar-slicing gauge. The physical situation was to describe the dynamical transition to spontaneous-scalarization and 
the further gravitational collapse to a black-hole formation. In Novak's works, the scalar-wave emission was analyzed 
and confronted with the sensibility of current detectors. 

Using a perturbative analysis, STT have been also analyzed in neutron-star oscillations and the emission and 
detection of scalar- waves have been confronted with the imprint left on the gravitational wave (spin-2) spectrum |l2j| . 

Except for the above studies (in spherical symmetry or using perturbation theory) up to the author's knowledge, 
there is no further attempt in STT for studying the dynamics of spacetimes with less symmetries or with another 
kind of matter (different from a perfect-fluid). 

Actually, one of the scopes for a future work is to analyze the dynamical transition to spontaneous scalarization 
and gravitational collapse with scalar- wave emission but for the case of boson stars (as opposed to neutron stars) . In 
this case the matter is represented by a complex scalar-field ip with the following energy-momentum tensor: 



i|VVf + W| 2 ) 



(94) 



TT = (V„V*)(V^) + (V6V*)(V o V0-fla* 

V^*i>) = my^ + A^) 2 , (95) 
with the corresponding 3+1 variables given as 

£matt = ^(|n V ,| 2 + |QV'| 2 )+'MIV'| 2 ) , (96) 



2 

nab r\ * a /-) b , /~)*br^a r^ab 



~(|Q*Mn*| 2 )+VM s 



(97) 



where 



5ma tt = i( 3 |n^-|Q^)-3^(M 2 ) , 

Smatt-^matt = |^ | 2 - | 2 - 4V^ (|^| 2 ) , (99) 

jr u = E^Ql* + %Q* . (100) 

Qt ■= Dai> , (101) 
V& := Cnip = n a W a ^ = A (d^ + N a Qt) , (102) 

and Q*^ and II* ^ are obtained by conjugation from the above variables. Here | | indicates the norm of the complex 
fields. 

From the point of view of the numerics, this is much more simpler than the realistic and difficult analysis of neutron 
stars. In this case, the problem involves an extra Klein-Gordon equation 

□V> = a^vv(^>) , (103) 

[which can be easily written in 3+1 form like Eq. I|71|l ] instead of the relativistic Euler equation for a perfect fluid. 
This problem involves therefore two scalar-fields (<j) and ip). 

In the Einstein frame and for the stationary situation, the phenomenon of spontaneous scalarization was shown 
to exist also in boson stars Nevertheless the dynamical analysis to that "phase transition" has not been yet 

performed. 
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VI. OUTLOOK AND CONCLUSIONS 



In this paper the field equations of scalar-tensor theories of gravity in the Jordan-frame representation have been 
recasted in a 3+1 form. This new system of equations shows that the Cauchy (initial value) problem is well formulated 
in that frame. This result opens the doors to several possibilities, both analytical and numerical. For instance, it 
would be interesting to write the 3+1 equations in full first-order form in order to perform a hyperbolic analysis along 
the lines of recent research in pure GR. In this direction, several gauge conditions for STT which generalize gauge 
conditions used in GR, have been proposed. These latter correspond to evolution equations for the lapse and the shift 
vector. 

Regarding the applications, a future plan is to use the formalism presented here to the case of several astrophysical 
phenomena. One of this is the numerical analysis of the transition to spontaneous scalarization and the collapse of 
boson stars. This would be similar but in many respects simpler than the analysis performed in neutron stars [lfj. 

On the other hand, when applied to specific spacetimes (for cosmology and spherically symmetric spacetimes) 
accompanied with perfect-fluid models, the new system of equations was shown to reduce to the equations analyzed 
numerically in the past by the author and colleagues [HI, Ha Ell • ^his illustrates that the 3+1 formulation derived 
here is successful in producing concrete results suitable for confrontation against observations. 

In this paper, it was also shown from the covariant form of the equations, that STT have a well posed Cauchy 
problem using a generalization of the harmonic coordinates (Appendix A). For completeness, the linear and the 
Newtonian limits of STT were re-analyzed in Appendix B. 
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APPENDIX A: ON THE WELL-POSEDNESS OF THE CAUCHY PROBLEM OF STT 

In the following, the proof of the well-posedness of the Cauchy problem of the system of Eqs. (|llf> and H17[) in 
the absence of matter will be considered. To do so, it will be convenient to chose a generalization of the harmonic 
coordinates x a given by the condition 

H a := Ux a - GV a (f> = , (Al) 

where Q is a function of the scalar-field which is to be determined so that the Einstein equations acquire a desired 
form. Here V a </> = g ab Vb4> and as in Sec. IV, is to be understood in the sense of four functions in this particular 
coordinate system. As before, 

nx ° = -La b [V^9 ab ] = d b g ab + l -9 ab 9 cd d b9cd = -g cd T a cd . (A2) 

Now, Eq. (|ll(l can be written in terms of T := T a a , as 

R ab = 8vrG (r ab - ^g ab T^j . (A3) 

Using the explicit expression for the Ricci tensor in terms of the metric and Eqs. IjAlfl and (|A2|> . one has 



Rab = R a i b -Qdi b 4>-g c{a d b) H c = MG [T ab -~g ab T\ , (A4) 

R% ■= -\g cd d 2 cd g a b + Kb{9,dg^,d4>) , (A5) 

where £f ab (g, dg, (f>, d(j)) are non linear functions of the metric, of the scalar field, and their first derivatives (indices 
are omitted for brevity in the arguments). 
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On the other hand the r.h.s of Eq. (|A3|) reads 

Rab = jdl^ + <% ab (g, dg, <j>, 8<f>) . (A6) 

where Eq. (|17|) was used to replace in terms of its r.h.s which contains no higher than first order derivatives of 
the scalar field, and also the relation T = S — E was employed (without the matter terms) which is given in general 
by Eq. $J$ . 

Therefore Eq. (|A3|) reduces to 

R% - 9c(ad b) H c =(g + i^j Sfat + ^ ab (g, dg, t/>, 90) . (A7) 
Q = -j , (A8) 



The convenient choice 



leads then to 



The gauge Eq. (|A1|I reads then 



Rub - 9c(ad b) H c = <% ab (g, dg, <j>, d<j>) . (A9) 



f 

H a = Ox a + jV a <j!) = Ox a + V a T = , (A10) 

where T = InF, which is exactly the same as the PHG 1)760 for = 1. 

Now, assuming that H c — can be maintained during the evolution (as shown below, this is indeed the case) , then 



Rab = m ah {g,dg,4>M) . (All) 

From Eqs. IjAlip and (|A5|) one has 

g cd dl d g ab = <€ ab {g,dg,<]>,dcj>) . (A12) 

On the other hand Eq. (|17|) together with the gauge Eq. I|A10|) . yields 

g cd d 2 cd ^ = ®{gA,d4>) . (A13) 

The non-linear functions £/ a b, 3§ab, ^ab, and & contains the explicit dependence of 4> from /, /', /" and V, V, and 
there appears also terms like 1//. Here it is assumed that the function / is a positive definite smooth function, and 
that V((f>) is also smooth and therefore that ff a b and 'Si are smooth functions of its variables. According to a theorem 
by Leray 0] a system of quasilinear partial differential equations like (|A12|) and i|A13fl has a well posed initial Cauchy 
value problem in the sense that given some initial data in E t and assuming (M, g ab ) to be globally hyperbolic, the 
system has a unique solution and the solution depends continuously on the initial data (see Theorem 10.1.3 of Wald 

What remains to be proved is that H a = in a neighborhood of St given H a = and dtH a = on E t . First, 
H a = on S f simply provides the initial values dtN and dtN 1 on E t given the set (hij, Kij, <j>, n)s t that satisfy 
the constraints n a R ab = 8irGon a (T ab — \g a bT) ■ Moreover, if the initial values are chosen as to satisfy the "reduced" 

constraints n a R^ b = n a &) ab (g, dg, (j), d(f>) on E t , this implies ( n a g c ( a d b )H c ) = 0, which in turn leads to (d t H a )s t = 0. 

In this way the pseudo-harmonic gauge (|A10)I and its time derivative are achieved initially. 



One recognizes in Eq. ElOl the gauge □*x a = where □» is the D'Alambertian with respect to the conformal metric g* b = F({f>)g a i,. 
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Finally, one can consider the effective Einstein tensor, G a b '■= G a b — F^ b TT , where G a b is the usual Einstein tensor 
and Tg TT is the energy-momentum tensor <|13|) in the absence of matter. The gravitational field equations can be 
written as 

G ab = 8wG eS T^ u , (A14) 

which in the absence of matter and in the PHG are equivalent to Eq. (|A11(I . Now, the key point in the proof for the 
PHG to be maintained during the evolution will be the use of the effective Bianchi identities V b (G a b/G c ff ) = 0. These 
latter result from the usual Bianchi identities and from the field equations. In fact, as mentioned in the introduction, 
these identities in the presence of matter lead to the conservation of the energy-momentum tensor of matter alone. 
However, in the absence of matter, which is the case at hand, the identities reduce simply to V h G ab = 0. One has 
then that the effective Einstein tensor can be written in arbitrary coordinates in terms of R^ b [which is given by Eq. 
(|A"5)) with ^ ab (g,dg,4>,d(/)) instead of gf ab (g,dg,(j),d(f>); cf. Eq. QA12J1 ] and H a given by Eq. (|A10j) as follows, 

G ab = Ri b - \R k 9ab ~ 9c(ad b) H c + \gabd c H c . (A15) 

Now, the use of V b G ab = and when Eqs. I|A12() and (|A13|) are satisfied, leads to 

V b G afc = -\g ab g cd dl d H h + S a {3H, ...) = , (A16) 

where $ a (dH, ...) includes lower order terms linear in dH. From the above equation one easily writes 

g cd d 2 cd H a = ,^ a (dH,...) , (A17) 

for which there exists a well-posed initial Cauchy value problem. Therefore, this equation ensures that H a = in the 
region where solutions of Eqs. (|A12|) and (|A13|I exist, if (^H a ^j = = (j) t H a ^j . In this way the pseudo-harmonic 
gauge is completely achieved. 



APPENDIX B: THE LINEAR LIMIT OF STT 



1. The linear limit "covariant" approach 

The linear limit of STT has been analyzed in the past by many authors (see Refs. 3]). For completeness, this limit 
is reanalyzed here and the point of departure will be the field equations in the Jordan frame. As usual, one considers 
first order perturbations of the Minkowski spacetime: 



9ab P 


- Vab ~ 




(Bl) 


Tab ? 


° 1 ab 




(B2) 




- 0o 4 


-ecf> , 


(B3) 


F{cj>) P 


^ F 4 


~eF$ , 


(B4) 


d+F{4>) P 




-eF^ , 


(B5) 








(B6) 




* Vol 




(B7) 






-eV^ . 


(B8) 



where t«l and the knot indicates quantities at zero order. In the 4+0 covariant formulation one can introduce the 
combination 

lab ■= lab + KVab'P , (B9) 
1 

lab ■= lab - -xVabl , (BIO) 
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where 7 = 7 a a and k is a gauge constant to be fixed later in order to simplify the equations. It is to be reminded that 
four dimensional indices are lowered or raised with the flat metric rj ab (or r] ab ). 
The resulting linearized Einstein Eq. is 

G ab = d c d(b^ a ) c - ^D^7afc - \^] a bd c d d lcd 



= 8ttG, 



of ab + k (dl b 4> - r) ab n n £j . (Bll) 



where T® b — results from the self-consistency of the perturbations at first order. Here is the D'Alambertian 
operator compatible with the flat metric rj ab . 

The following gauge condition will be imposed in order to simplify the linearized equations: 

d c % b = . (B12) 



As it is shown below, for the choice of k given by Eq. (|B20|) , this gauge condition corresponds to the linear limit of 
the pseudo- harmonic gauge Eq. I|76(l for 6 = 1. 

Then, from Eq. i|Bll|) the resulting wave equation is 



□„7«i = -16ttGoT o6 - 2k (d 2 ah 4> - rjabO^ . (B13) 
The linear approximation of the effective energy-momentum tensor Eq. l|13f) and the Eq. l|17f) turn to be 



T^matt f/ 
rfi ± ab 1 JO 

In.; 



ft*™. ■ ^rS^-"^'*) • (B14) 

n v $ - mU = 47ra^f mat t , (B15) 
Jo 

<■■=— %rr , (B16) 

/ 3f' 2 \ ' ( B1? ) 

where the following conditions were used 

V{<h) = = V ' = , (B18) 

resulting from the consistency at first order of the linearized Einstein and scalar-field equations and assuming that 
^matt ab — 0- ^ n ^is way the wave equation l|B13|> becomes 

2 ~ f ft, 

T^matt of,.. 1 



□„7«& = -jTT 1 Vab^ ) . (B19) 



Note that at this order of approximations the application of the ordinary divergence in Eq. (|B19|) and the use of the 
gauge Eq. (|B12|) leads to the energy-conservation of the matter perturbations: d a T£ att = 0. 
The convenient choice for k, 



Ik 
fo ' 



(B20) 

allows to simplify Eq. I|B19(1 considerably. In summary, one deals with the following wave equations for the gravita- 
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tional and scalar modes 10 : 

^u b = -T f T U , (B21) 
Jo 

d a j ab = , (B22) 



- m 2 4> = 47ra^f matt , (B23) 
Jo 



with the constants TOq,q; and k given by Eqs. (|B16(I . I|B17|I and (|B20(I respectively. The analysis of propagation of 
gravitational and scalar waves will be not pursued here, this has been done elsewhere (see Refs. Q ). 

2. Linear limit of 3+1 equations 

In order to linearize the 3+1 set of equations (5$, {2jJ), and ©-JZU, 

once again assume first order deviations 

of the Lorentz metric r) a b as follows 

N w 1 + dV , (B24) 
iVi « eiV, , (B25) 
hj « 5y + e/iy , (B26) 

where e< 1. 

Note that at first order the inverse of the 3-metric is given by 

« % - ehij , (B27) 

such that the condition huh 1 ^ = 5 l j holds. Therefore, covariant and contravariant components of tensorial quantities 
having no zero order terms are identical to each other at first order. For instance N l — h ll Ni w eJVj. In order to 
compare with the 4+0 (covariant) linear approximation 



it turns out that 



Qab ~ Vab + ejab , (B28) 

700 = -2N , (B29) 

7o* = -Ni , (B30) 

7ij = ^if • (B31) 



The 3-Christoffcl symbols turn to be 



z T) k « e 3 f} fc , (B32) 

3 f jfe := ^(-flihifc + flfcfttf+fliftw) • (B33) 

Therefore 3-covariant derivatives of 3-tensors having no zero order terms become at first order ordinary derivatives. 
For instance 3 V-,A r i w edjNi. Then Eq. ggj leads to 

iCij « ei^ , (B34) 
:= —(dthij+WiiNfl) , (B35) 



It is not a surprise to recognize in -fat, the conformal metric perturbation 7*j/i*b which is defined from g* b ss Fo + e ^7 a 5 + -p^-^ab 
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with 



K' 



eK 



K := K L 



eK , 



K := -1 (d t h + 2diN l 



(B36) 
(B37) 

(B38) 



where h :— ha is the trace of the 3-metric perturbation. Using the definition of the 3-Riemann tensor in terms of 
the 3-Christoffel symbols it is easy to obtain the linearized approximation for the 3-Ricci tensor and the 3-curvature 
respectively: 



i R l 



e 3 R 



V i 
2 



e 3 R , 



'R 



3 R := d 2 kl h kl - 3 Ah 



(B39) 

(B40) 

(B41) 
(B42) 



where 3 A := d\ stands for the Euclidean 3-Laplacian. 

Concerning the scalar field and the matter, in addition to the approximations given by (|B3|t — (|B8|) . the following is 
also to be assumed 

(B43) 
(B44) 
(B45) 
(B46) 
(B47) 
(B48) 
(B49) 

where the knot indicates quantities at zero order. 

The Hamiltonian constraint Eq. (|67|l when linearized reads 









Qi 


:= di<f> , 
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« en , 
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:= d t 4> , 




^ matt 


TTimatt 1 
W ii 4 


^^imatt 


Tinatt 
J i 


jO matt 


+ e att 


omatt 
D ij 


^ qO matt 


+ eS™ att 



2E n 



or explicitly 



fo 



d 2 kl h kl - 3 Ah 



fo 



Vo 
fo 



diQ 1 



2E„ 



fo 



The linear limit of the momentum constraint Eq. H68|) is 



diK 1 



d l K + fy-difl 
Jo 



J? 



Finally the dynamic Einstein Eqs. (|69|) when linearized read 

fo ; 



d t k\ ,+ &■ a-iv - 3 r i 



foV "4> - (S n 



fo 



d l Q 3 



2/o(l + %- 



fi 2 
Jo 

fo 



(B50) 



(B51) 



(B52) 



2S, 



matt j 



Mi 
2/0 



(B53) 
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The linear limit of the evolution Eq. I|7(JII reads 



2 



d t K + *AN - fftQ 1 = _■ i SUtt + ^matt 1 + - • (B54) 

/o 2/0(1 + ^) I V /o / J 

Finally, the linear limit of the Eqs. (|31() and l|71|) takes respectively the following form 

d t Qi = din , (B55) 



o tt n Ai /o^ + 2/o^inatt 

9 t II - djQ 3 = -ppr — . (B56) 



Equations l|B44|l . (|B46|) and (|B56|) are indeed equivalent to Eq. (|B23f) . One reminds T mat t = .S^att — E maM . 

It is to be emphasized that in the above linearization the self-consistency of the 3+1 equations up to first order 
imply that the zero order source fields must vanish identically: 

^matt = = jO matt = gC I matt = yt _ y Q (B57 ) 



The above equations (|B51fl — (|B53|) are the 3+1 decomposition of the 4+0 equations (|B11|) 
It is to note that the 3+1 splitting of the perturbations (|B9|) is given by 



700 = i(£-2/c£-2JV-) , (B58) 

7o* = -Ni , (B59) 

% = hij - (h - 2n4> + 2JV) , (B60) 

7 = - (h + 27V) + 4k0 . (B61) 

Now, assuming the PHG given by Eqs. (|9"2|) and (|9"5|l and taking n — —Q!F' Q = — 0/o//o, the linear limit reads 



d t N + kIL + K = , (B62) 
1 

2' 



d t Ni + djhij — diN — -dih + nQi = , (B63) 



where the following linear approximations were used V 1 = -^dj (y/hh^J ~ eVi with Vi = — [djhij — ^dihj. This 

result is obtained using Eq. (|B27|) and also from the fact that at first order h = de thjj » 1 + eh. Both Eqs. (|B62|) 
and (|B63|> amounts to the 3+1 decomposition of the Lorentz-Dc Donder gauge Eq. (|B12|) . 
From Eqs. (|B62|) . (|B54fl . and (|B56|) one obtains a wave equation for N: 



(B64) 



n v N -(k + &) diQ 1 = 1 { 5 matt (1 - K f ) + E matt [1+VL + K A - 2f V^ (k + §^ 

v fo j ■>/,,( \ ■ %-\ 1 v ■ / » / v 2./, 

where O v := -dj t + 3 A. 

On the other hand, the combination of Eqs. HB52() and (|B35|> with Eqs. (|B63|) and (|B62|) (noting d t Q l — which 
follows from the integrability condition a^cj) = df t <f>), results in a wave equation for N l : 

U v N l ~2 (n+ff) d t U = -^f^ . (B65) 

\ JO J JO 
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Finally, the linearized dynamic Einstein Eqs. (|B53|) together with Eqs. i|B35(l and 1B63(I lead to a wave equation for 



hij as follows, 



+ ( n+ I d t Q 3 



fo 



2/o 1 



f,v"<>- i^n.n -■ £■„„.,.) ( i - + 4-J 



25" 



3/r 

2/0 



From Eq. (|B66|I and tracing one obtains 

, ft 



fo 



diQ 1 



-Smatt + 3/oVg'V + 3-E ma tt 1 + 



fl 2 

fo 



(B66) 



(B67) 



Therefore combining Eqs. (|B64|) and (|B67|) yields 



fo 



□„ [ ~-N) + 2[K+^)d l Q l = - 



2/o(l + ^ 



S, 



matt ^matt 



t ) - 2 K / f " 



4£ n 



3/r 
2/0 



Finally one can further combine the wave equation l|B56[) and (|B68|) to find 



N — K<j) I + 2 ( K + 



/o 



2 -Err 



/o 



(B68) 
(B69) 



Perhaps a more straightforward way of obtaining this wave equation is by applying dt and d l to Eqs. I|B62() and (|B63|) 
respectively and then adding the resulting equations followed by the use of Eqs. (|B51I) and HB38I) . 
On the other hand, the combination as given by Eq. (|B60(I provides the wave equation 



n v 7 y - + 2 [k 



fo 
fo 



2 Cjrnatt 

bQ, = f— 

Jo 



fo I •S'matt — Ematt ) + 2/oVq"0 



(B70) 



Therefore Eqs. (|B65|) . I|B69|I . and i|B70(l are equivalent to Eq. <|B19(1 . In such equations one can employ the choice 
k = — /0//0 (which corresponds to 9 = 1) to simplify the expressions as follows 





2-Ematt 
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(B71) 
(B72) 
(B73) 



These equations are equivalent to Eq. i|B21(l . 

In fact, with k = —fo/fo, another useful wave equation is obtained from Eqs. I|B64(1 and (|B23|) . 



□„ N 



2/V 



matt "T -Ematt 



2/0 



(B74) 



3. The Newtonian approximation 



Since the complete equivalence between the linear limits of the covariant field equations and the 3+1 equations has 
been shown, the Newtonian approximation (slow varying sources and small pressures) will be considered only from 
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the linear limit of the 3+1 approach. In order to do so, one takes -E ma tt = P , S matt <C p, J matt <C pc and neglect the 
time dependence of the perturbations. From Eqs. (|B72|) and (|B73|) where one assumes 5™ att « « ./matt an d taking 
into account that perturbations are well behaved at infinity, one concludes that the unique solutions are N l = const, 
and 7jj = const. Such constants can be gauged out and so iV l = = 7^ . 
On the other hand, Eqs. I|B23(1 and l|B74(l lead respectively to 
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(B76) 



The Newtonian approximation in the geodesic equation for test particles leads to the identification of N with the 
Newtonian potential <&at. Hence, from Eqs. (|B75|I and l|B76|) one finds the solution 
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8^/0 



a (ft 
2 \fo 



p(x' ) e ~ m ° lS ~ S ' 1 d 3 x , 



B.C. 



(B77) 



where to and a are given by Eqs. (|B16|) and (|B17|) . respectively. 

For instance, in the exterior of an extended body with mass-density p with compact support, the Newtonian 
potential is given by 
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(B78) 



where r = \x\, and /o = Fq / (8ttGq) was restored. In the above expression M corresponds to the gravitational mass 
of the body given by the volume integral of the mass-density p . 

It is customary to express Fq in terms of the effective Brans-Dicke parameter Eq. (0J 
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(B81) 



Different experiments bound the parameters 
tonian force law ~ 1/r 2 (see for instance Refs. 



) of theories that induce modifications of the New- 
)■ 



Now, for the massless case to,q = 0, the Newtonian potential reads 
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(B82) 



where 
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(B83) 



is the dressed gravitational constant measured by a Cavendish experiment 0, 0] . 



In the strong-field regime (e.g., neutron stars) and for instance in the static and spherically symmetric case, the asymptotic behavior 
for the lapse and the scalar-field are N ~ 1 — GoMadmA' and 4> ~ uifi/r (for mo = and <f>o = 0), where is the scalar charge llll 
and A/adm is the ADM mass. 
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APPENDIX C: 3+1 DECOMPOSITION OF DIFFERENTIAL OPERATORS 

Let us consider the D'Alambertian of any number A of functions ip (x a ) with their corresponding source terms: 

u^ A = S^ A . (CI) 

Now, 

¥ A = ~^=dc [V~99 cd d d ^ A ] . (C2) 



g 

Using g ab = h ab — n a n b and the fact that g := detg a b = —N 2 h where h := det/iy, then 

□V> A = —^dc \NVhh cd d d ^j A ] l —d c \NVhn c n d d d ip A ] . (C3) 

Hence, 

Di/f A = 3 Ai) A + a c V c ij) A + Kn c V c i/) A - n c V c (n d \7 d ip A ) . (C4) 

where one reminds 3 A is the Laplacian compatible with the 3-metric hij and a c — n d V d n c = D c [\nN] is the 4- 
acceleration of the observer with 4- velocity n a . Here it was used the identity K = — V c n c . In order to obtain a system 
of first order equations one can further define as in Eqs. (|26|) and l|27|) the following quantities 

Q aA : = D a ip A , (C5) 
IT 4 := £ n </> A = n'VcT// 4 . (C6) 

Collecting the above results one obtains 

£ n Il A - a c Q cA - D c Q cA -T1 A K = -S^ A . (C7) 
A simple application of the above results is the case when Tp A = x A (with x A — (t 7 x l ) as in Sec. IV). It turns then 
IP 4 = n A = ±{l,N l ) and Q A = h d 5 A = h A [where it is to note that Q c ° = and D c Q ci = -^d 3 (Vhh^ = V% 

The above Eq. I|C7[) with a suitable choice of S^ A is then equivalent to Eq. (|82|l which provides an evolution equation 
for the lapse and for the shift. For instance, taking S^ A = in Eq. I|C7|) allows one to recover Eqs. (|92() and (|93|l (for 
= 0) which are equivalent to the harmonic coordinate condition Eq. 175(1 . 

On the other hand, one can take U^ A = with (no sum over index A). 

Now, using n c n d — h cd — g cd , one obtains 

n c n d W c W d tp A = h cd V c V d ilj A - U^ A . (C8) 

Using the orthogonal decomposition (see Sec. Ill) V a -0 A = D a ip A — n a n c '\7 c ip A — Q A — n a TL A it yields 

n c n d W c W d ^ A = D c Q cA + KH A - D^j A , (C9) 

where one used /i cd V c n d = V c n d = ~K and h cd n d = 0. In this way the equation Otp A = becomes 

^ip A — ■ * [TC> (D c Q cA + if IP 4 ) - GV A F] . (CIO) 
L + L^a 

Finally, — d^p A is given by the l.h.s of Eq. I|C7(I which leads to 

C n T\ A - a c Q cA - D c Q cA - U A K = — 1 [Tf> (D c Q cA + KIL A ) - QV A T] , (Cll) 
which simplifies to give 

C n Tl A - a c Q cA = - 1 (D c Q cA + KT1 A + QV A T) . (C12) 

1 + T.^A 

As before, with tJj a = x A = (t, x l ) and T^a = T qA , and using Eqs. (gSJ and |2HJ (which leads to Ii A = n A = ± (1, N l ) 
and Q A = h d 8 A = h A ), from the above equation (|C12|) one recovers Eq. (|82|l for B — 0, A = 1. Following the above 
examples and taking for S^ A the r.h.s of Eq. I|81|) one can then generalize the expression of Eq. IJC12I) so as to obtain 
Eqs. (H3> and (JB3J. 
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APPENDIX D: SELF-CONSISTENCY OF 3+1 FIELD EQUATIONS 

As mentioned in Sec. I, one of the main properties of the STT is that they preserved the weak equivalence principle 
in the sense that the conservation of the effective energy-momentum tensor (|13|) 

V b T ab = , (Dl) 

together with the field equations and 112|) lead to the conservation of the energy-momentum tensor of matter 
alone 

V 6 T« b att = . (D2) 

The same result can be expressed in terms of the 3+1 decomposition of the field equations as performed in Sec. III. 
However, it is not at all evident that the same result will turn out if one has committed any mistake (one as "innocent" 
as an error sign) during the process of the derivation of the 3+1 field equations of STT. In order to carry out such 
a check one requires to express the conservation equation (|Dlll in 3+1 form. Without showing any mathematical 
details, the projection of Eq. IjDljl on T, t and using Eq. leads to the equation of conservation of momentum 

h a c n b W b J c + D b S ba - (K a c + Kh a c ) J c + (S a b + Eh a b ) D b [lnN] = . (D3) 

Of course this can be written as an evolution equation for J c by noting h a c n b V b J c = h a c CnJ c + K a c J c . 
In a similar fashion, the projection of Eq. (|D1|) on n a , leads to the equation of conservation of energy 

n b V b E + D b J b + 2J b D b [\nN]- S ab K ab - EK = , (D4) 

which provides the evolution for E since n b W b E = LnE. 

The next step consist in using the expressions (|62|) — <|64|) in the equations (|D3(1 and i|D4|) . It is clear that the 3+1 
field equations will be required in order to simplify the equations. An extremely involved calculation that demands 
the use of the evolutions Eqs. J22J, (|3*T)l . (|rJ$|l and lj7T)l as well as the constraints Eqs. (|u7|) and ljrl%)l leads finally to 
the desired result 

h a c n b V b J c mM + D b S b : att - (K a c + Kh\) J£ att + (5» att b + E matt h\) D b [\nN] = , (D5) 



n b V b E matt + D b J b natt + 2 J 6 a « A, [InA] - S± tt K ab - E matt K = . (D6) 
This self-consistency check ensures that the 3+1 field equations of STT in Sec. Ill are correct in their full form. 
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